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Abstract: We examine the phase structure of massive Yang-Mills theory in 1+1 dimen¬ 
sions. This theory is equivalent to a gauged principal chiral sigma model. It has been 
previously shown that the gauged theory has only a confined phase, and no Higgs phase in 
the continuum, and at infinite volume. There are no massive gluons, but only hadron-like 
bound states of sigma-model particles. The reason is that the gluon mass diverges, being 
proportional to the two-point correlation function of the renormalized held of the sigma 
model at x = 0. We use exact large-IV results to show that after introducing a lattice reg¬ 
ularization and typical values of the coupling constants used in Monte Carlo simulations, 
the gluon mass becomes finite, and even sometimes small. A smooth crossover into a Higgs 
phase can then appear. For small volumes and large N, we find an analytic expression 
for the gluon mass, which depends on the coupling constants and the volume. We argue 
that this Higgs phase is qualitatively similar to the one observed in lattice computations at 
N = 2. 
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1 Introduction 

Yang-Mills theories coupled to a matrix-valued Higgs field are known to have a confined 
phase and a Higgs phase in 2+1 and 3+1 dimensions [1], The phase diagram in 1 + 1 
dimensions is much less understood. 

The author and P. Orland proposed studying the (l+l)-dimensional theory as a gauged 
principal chiral sigma model (PCSM) [2], We found that in the continuum theory at 
infinite volume, there is only a confined phase. The physical excitations are hadron-like 
bound states. We computed analytically the spectrum of the meson-like hadrons in the 
nonrelativistic limit, using knowledge of the exact S-matrix of the PCSM. 

A numerical lattice study of the phase diagram of the same theory, with gauge group 
SU(2), was published by S. Gongyo and D. Zwanziger [3]. In this reference, they found 
evidence for both the confined and the Higgs phase, even in 1 + 1 dimensions. There seems 
to be a crossover between the two phases, instead of a sharp phase transition. This appears 
at first sight contradictory to the results of [2]; however, what they found is that the 
crossover seems to disappear as the volume of the system is increased. Their calculations 
suggest that there is a Higgs phase, but it disappears as they move towards infinite volume. 

In fact, within the paper, the authors make the statement, “from this data at finite volume 
we cannot conclude that there is a phase transition to a symmetry-breaking [Higgs] phase 
at infinite volume". Furthermore, as we will suggest, the value of gluon mass monotonically 
grows as the continuum limit is taken. 

Inspired by the nontrivial results of [3], we extend the approach started in [2] to study 
analytically the (l+l)-dimensional theory with a lattice discretization and at finite volume. 
Our analytic computation is so far only possible at large-IV, so it cannot be directly com¬ 
pared to the N = 2 results of [3]. Nevertheless, we will see that the behavior of the Higgs 
phase we find at large N is qualitatively similar to the N = 2 results of [3]. 

The mathematical object we study is the two-point correlation function of the renormalized- 
field operator of the PCSM. This function has been calculated at infinite volume in [4], [5], 
and finite volume in [6], using the LeClair-Mussardo formula for integrable field theories 
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Our main claim is that even though there is only a confined phase in the continuum 
theory at infinite volume, a Higgs phase can arise once a lattice regularization is introduced, 
particularly at very small volumes. Our calculation is valid only for large N, since the 
nonperturbative results needed are only available in this limit. We will later argue that 
a similar phenomenon can occur at finite N (particularly at N = 2, which is desired for 
comparison with numerical results), however, we can only rely on perturbative calculations, 
and only study the infinite-volume limit. 

The PCSM has the action 

'S'pcsm = [ d 2 x 7 ^Trd fl U\x)d fl U(x), (1.1) 

J ^9 o 

where U(x ) G SU(1V). This model has been shown to be integrable, and its exact S- 
matrix is known [8]. The action (1.1) has an SU(./V) x SU(A') global symmetry given by 
U(x) —» VlU{x)Vr , with Vl,r £ SU(1V). The PCSM is asymptotically free and has a mass 
gap, which we call m [9]. 

One can obtain the massive Yang-Mills action by promoting one of the SU(fV) sym¬ 
metries of (1.1) to a local gauge symmetry. We choose to gauge the left handed symmetry, 
Vl —> Vl(x). We introduce the covariant derivative D^ — ieA^^x), and field strength 

F /tu = d^A v — dyA^ — ie[A^,A u \. The gauge held, A^x) transforms under the gauge 
symmetry as A M (x) -» V^ l {x)A^{x)Vl{x) — ^vl{x)d^VL{x). The gauged PCSM action is 


S 


d 2 x 


--TAF^F^ + -^Tr (D^D^U 

4 z 9 o 


( 1 . 2 ) 


In Reference [2], it was argued that the model with action (1.2) is always in a confined 
phase, rather than a Higgs phase. The argument was made by examining this action in the 
axial gauge A\ = 0. In this gauge, the action (1.2) becomes 


S 


d~x 


-Tr^Mo) 2 + -^Tr(5 

^ ^9 o 


-^Trdit/^i U 


The gauge held Aq can now be integrated out, obtaining the action 


s = I d2x • (L3) 

where i^(®)fe = —iTi'ti > d fi U(x)U^(x) is the Noether current associated with the left handed 
symmetry, and tf,, with b = 1,... ,N 2 — 1, are the generators of SU(1V). The action (1.3) 
describes a PCSM with an additional potential between the left-handed charges. 

Looking at (1.3) one would naively conclude that, since U'U = 1 (because U is unitary), 
the potential between the charges is screened by a force carrying excitation of mass e/g o- 
This would be interpreted as the massive gluon of the Higgs phase. As was discussed in 
Ref. [2], this reasoning is wrong. The main problem is that even though the bare held 
U{x) is unitary, the physical renormalized held, which we call d>(x), is not [4], The relation 
between the two helds is given by 


(0|Tr$(x)T(0) t |0) = Z[ 5o (A),A]- 1 (0|TrC/(x)I/(0) t |0), 
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Where Z\go(A), A] is a renormalization constant, and A is the momentum cutoff of the 
theory. Because of the asymptotic freedom of the PCSM, the renormalization constant 
vanishes logarithmically with the cutoff. This means that renormalization pushes the mass 
of the gluon to infinity. Therefore one does not see any massive gluons in the physical 
spectrum of the theory. The model is actually in a confined phase. 

Using the fact that for the bare field, t/t(0)C/(0) = 1, the renormalized gluon mass, A4, 
is given by 


M 


-Zbo(A),A ]- 1/2 

9 o 


- [^(0|Tr<h(0)<f>(0)t|0) 
9o [A 


(1.5) 


The only reason that there is no Higgs phase is that the PCSM field-renormalization con¬ 
stant vanishes as we increase the cutoff. Consequently, the gluon mass diverges. It is then 
clear that once a cut-off is introduced by placing the theory on a lattice, the gluon mass 
becomes finite. In this letter we propose that placing the model in a small volume can 
further suppress the divergence. Given the values of the couplings, e, go used in [3], the 
mass becomes small enough to be measured in Montecarlo simulations. We show that the 
dependence of the gluon mass on the volume agrees qualitatively with the results of [3]. 


2 Infinite volume PCSM correlation function 

An expression for the two-point function of the renomalized field of the PCSM was found in 
Ref. [4]. This expression was found using the fact that the PCSM is integrable, combined 
with the large-IV limit. The approach is to first find the exact S-matrix, and all the form 
factors (matrix elements of local operators) of the renormalized field using the integrable 
bootstrap program [10]. Once all these form factors are known (they are known for the 
PCSM only in the large-IV limit, with g^N fixed), the two-point function is given by the 
spectral sum, 

M®) = jf E^°l $ ( X ')^o[^oao(0r|0) 

aofio 

= v £ £ (oi*(°)*.«.i») <*i i#(o)*»«r io>, (2.D 

ao,bo 'I 1 

where |T) is any state with particles and antiparticles, and pq, is the sum of the momenta 
of the excitations of the state IT). 

After introducing the exact form factors in (2.1), the exact expression for the two-point 
function given in [4] is 

w(x) = L Mi -L <i " 2wexp - •") n )*+**• (2 - 2) 

where the integration variables, 0j , correspond to the rapidities of the intermediate particles 
and antiparticles. These rapidities parametrize the energy and momentum of a particle by 
Ej = m cosh 9j, Pj = m sinh 0j. 
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For the purposes of this letter, we are interested only in the divergence at x = 0 of 
(2.2). If we set x = 0, the expression (2.2) diverges because one has to integrate over all 
the physical values of the rapidities. This divergence can be regularized by introducing a 
rapidity cutoff, —A < 9j < A. The rapidity cutoff, A, is related to a standard Euclidean 
momentum cutoff, A, by 


A = sinh 


-l 



In 


A / ^ - - + A / ^ + - \ 

V 2 m 2 2 V 2 m 2 2 J 



The regularized two-point function at x = 0 is 


i 00 pX pX 21 

= ih... n 


1=0 


L I (9j - 63 + 1) 2 + 71-2 


(2.3) 


The short-distance behavior of this PCSM two point function was first studied in [5], 
and the regularized function at x = 0 was shown in [6]. We will not show the full compu¬ 
tation, but only quote the result we need, and refer the reader to the original paper. The 
result of [5] and [6] is (after redefining integration variables Uj = 9j /A) 


W A (0) 


8n 2 



dmp n (u) 


2 

a n l + 


C 2 X 2 + 0(X). 


(2.4) 


where a n and <p n (u) are the eigenvalues and eigenfunctions of the fractional Laplacian oper¬ 
ator, A 1 / 2 = \J—d 2 / du 2 . The upper bound C 2 < 0.0219 was found in [5] 1 . The correlation 
function diverges, being proportional to (In A) 2 , which agrees with the perturbative result 
and asymptotic freedom [11]. 


3 Finite Volume Correlation Function 

The PCSM two-point function at finite volume has been calculated in [6]. This is done 
by compactifying the x 1 dimension into a circle of length V. In 1+1 dimensions, this is 
equivalent to placing the system in a finite temperature T = 1/V (if we compactify x° 
instead of x 1 ). Using the exact S-matrix, and the thermodynamic Bethe ansatz (TBA) [12] 
one can calculate the partition function. 

It was shown in [6] that the the partition function of the PCSM coming from the TBA 
at large N is trivial and equivalent to a free bosonic gas. However, the correlation functions 
of operators are not those of the free theory. 

The two-point function of the renormalized held was evaluated at finite volume using 
the Leclair-Mussardo formula [7], which we will not discuss here. This is essentially a 
spectral sum, similar to (2.1), but where one sums not only over all the intermediate states 
of the two-point function, but one also includes an ensemble average, with a second sum 
over states weighted by the thermal distribution function. 

1 We have learned by private communication that Eytan Katzav and Peter Orland have calculated C 2 

and found it is exactly y(—. We will only use the upper bound found in [5], as the calculation of the exact 
value has not yet been published. 
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There are several problems with the Leclair-Mussardo formula for two-point functions 
that make it not generally valid for all integrable theories [13]. The strongest objection is 
that generally, the expansion is not well defined, as there are poles for real values of the 
rapidities that need to be integrated over. However, it was argued in [6] that these problems 
do not affect the PCSM at large N, and therefore the two-point function should be valid. 

For this model, the poles do not lie in the real line of rapidities, and the Leclair-Mussardo 
formula is well defined. 

We are interested in the two-point function at x = 0, but in a finite volume, V. This 
result is found in Eq. (6.6) of Ref. [6]: 

.. oo 2 1 21—n 1 2l—n 1 —n 2 

wW = s EE E E dd i---/ d0 2/+i[/i(0i)/i(^+i)]^ 

1=0 n 1 =0 n 2 =0 n 4 =0 ^ ^ 

n 1 n^+n 2 n 1 +n 2 +n 4 21 

xn^i(M J+ i) n F 1 ,. 1 (d j , 0 j+1 ) n F- lt - 1 {e j ,o j+1 ) n 

j =1 j=n 1 +1 j=n 1 +n 2 +1 j=n 1 +n 2 +n A +1 

. oo 21 21—n 1 21—n 1 — n 2 

+ S EE E E dd!... d9 2l+1 [f. l {9 l )U{9 2l+l )}^ 

1=0 n 1 = 0 n 2 =0 n 3 = 0 ^ ^ 

n 3 n 1 +n 3 n 1 -\-n 2 -\-n 3 21 

x] jF-i >1 (e j ,e j+1 ) n F ltl (e j ,d j+1 ) F li . 1 (e j ,e j+1 ) J] F. h _ 1 (9 j ,9 j+1 ). 

j —1 j=n 3 -\-l j=n 1 -\-n 3 -\-l j=n 1 -\-n 2 +n 3 +1 

(3.1) 


where F au<Jj {9i,9j) 


{e. l -e j ) 2 +2 Wi ~ a ^-K'' 


with cTij 


±1, and f±i(9j) 


\J[\-\- gT^mcosh 6j\ 


The expression (3.1) is analyzed for both large and small values of V in [6]. At large 
V. one simply recovers the infinite volume expression (2.3). The small-volume limit is 
studied by realizing that for small volumes (for mV « 1), the functions ,/'±i (Oj) become 
approximately 


fm 


' i, -£<<?<£, 

1, otherwise, 


/- i( 0 ) 


' -C<9 <£, 

0, otherwise. 


where £ = In • 

Using techniques similar to those of [5], it was found [6] that for extremely small 
volumes, with C ~ A, 


W x {0) v « C 2 ( A 2 - C 2 ) + -C + O(C 0 ) + 0{ A - C). 

7T 


(3.2) 


At the smallest possible physical volume of V = 1/A the (In A) 2 term completely cancels 
out. For these very small volumes, the divergence of the two-point function is reduced from 
order A 2 to order A. 
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4 The Higgs phase and gluon mass 


We have presented the two-point correlation function of the renormalized field operator of 
the PCSM at large N in a finite volume. For very large volumes, the standard, infinite- 
volume two-point function from [4],[5], which diverges as (InA ) 2 at short distances, is 
recovered. For very small volumes, we saw that the order of this divergence reduces from 
(In A ) 2 to In A as we reduce the volume size. 

We would like to discuss in more detail what we mean by very small volumes. Suppose 
the momentum cutoff originates from placing the theory on a lattice with spacing a 1 /A. 
Now suppose the length of the x 1 direction is V = na. For the lattice computations of 
Ref. [3], for example, the values used were n = 256, 512, 1024. As an arbitrary example, 
for the sake of illustrating our point, let’s imagine a lattice spacing is chosen such that 
ma = 10' 14 (this is the assumption that m is small enough that all these volumes are still 
in the regime where mV < 1). It will soon be clear why this particular value of ma is 
chosen and is relevant. For this spacing, the rapidity cutoff is A = 32.236. We can calculate 
£ for the different volumes n = 256, 512, 1024, for which we find £256 = 26.691, £512 = 
25.998, £1024 = 25.305. That is, for typical lattice volumes, like those of Ref. [3], £ and A 
are of the same order of magnitude. The “softening" of the divergence from (3.2) is a large 
and significant effect for typical lattice computations at small volume. 

The question remains, is the PCSM mass gap, m, small enough that mV « 1 for the 
volumes studied in Ref. [3] ? The only length scale in the PCSM is m, and this parameter 
is usually fixed, and the bare coupling go (A) varies as we change the cutoff, A. Their 
relationship is [ 11 ] 


r ^ = ma = — (e ~ 2 */ 9 0 (4.1) 

A g 0 V / 

where K is some non-universal constant, which depends on the regularization procedure. 
The parameter that was controlled by hand in Ref. [3] is 7 = 1 /gf } . The values used were in 

the range 7 G {1.10}, with a crossover into the Higgs phase detected around 7 = 5 , 6 , 7, 8 . 

As is seen from Eq. (4.1), the mass gap is very small for the values of 7 at which the crossover 
was seen (and thus the size of the PCSM particles is very large). For example, if we take 
K = 1, we find 7710(7 = 5) = 5.07834 x 10~ 14 , and 7770(7 = 8 ) = 4.18334 x 10~ 22 . It is then 
clear that the computations of [3] are all done deep in the small-volume regime. 

We stated that the only reason for the absence of a Higgs phase in the continuum 
theory is that the gluon mass is proportional to the renormalized-field two-point function 
at x = 0. The gluon mass becomes finite once a momentum cutoff is introduced by placing 
the theory on a lattice. Depending on the value of the coupling constants and the volume, 
this mass can become small enough to be observed in Monte Carlo simulations. 

We can estimate the gluon mass using Eq. (1.5) and (3.2), with C 2 ~ 0.0219. We 
compare our mass estimates with those of Table 2. of Ref. [3]. The gluon mass is a 
function of the variables, go, e and n = V/a. We now switch to the notation of [3], where 
the couplings used are 7 , and [3 = 1/e 2 . For small volumes, the gluon mass is approximately 
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(with K = 1) 


M(j, /3, n) 






7 8 

T 



.(4.2) 


Equation (4.2) is not expected to agree quantitatively with the results of [3], as it is valid 
only in the large-lV limit. However, we argue that it qualitatively explains some of the 
observed behavior. 

First we examine the /3 dependence. As can be observed in Table 2 of [3], increasing 
the coupling, /3 while keeping 7 and n constant, generally decreases the gluon mass (this 
becomes clearer in their computations at higher values of 7 , closer to the continuum limit). 
The clearest example shown in [3] is for 7 = 8 , n = 256, where they find = 1-288. 

Plugging these same parameters into (4.2), we find = 1-291) which we believe 

confirms the /3 dependence of the gluon mass. 

The size of the lattice spacing relative to the PCSM mass gap is controlled solely by 
the 7 parameter (Eq. (4.1)). The continuum limit is equivalent to the limit, 7 —> 00 . It 
is clear from (4.2) that the gluon mass diverges as expected in this limit. For high values 
of 7 , keeping (3 and n fixed, we expect the gluon mass to increase as we increase 7 . This 
is confirmed in [3], where the clearest example is for their largest measured values 7 , fixing 
f3 = 120, n = 256, finding = 0.928. From the formula (4.2), for these values, we 

hud = 0-714401. Here we expect that the large -N limit plays a larger role, since 

the powers of the logarithms in (4.2) will be different at smaller N. 

Finally, we verify qualitatively the volume dependence of the gluon mass. As can be 
observed from Table 2 of [3], how the gluon mass reacts to a change of volume depends 
on the value of 7 . For small values of 7 , increasing the volume, slightly decreases the 
mass. For example, fixing /3 = 120, 7 = 5, they found A i(n = 256) = 0.240(2), and 
A4(n = 512) = 0.237(3) (in lattice units). However there is a change in behavior somewhere 
between 7 = 5 and 7 = 6 . For large values of 7 , increasing the volume slightly increases 
the mass. For example, for (3 = 120, 7 = 6 , they find A4(n = 256) = 0.265(5), and 
A4(n = 512) = 0.267(2). Our formula, (4.2) reproduces exactly this behavior. For a fixed 
value of /3 and 7 , the mass is modified very slightly by changing the volume from n = 256 to 
n = 512. For small 7 , the mass decreases with increasing volume, and for large 7 , the mass 
increases with volume. This change in behavior for Eq. (4.2) happens around 7 ~ 10.38. 
To illustrate with some examples, for /3 = 120, and 7 = 5, we have A4 (n = 256) = 1.71503, 
M(n = 512) = 1.70266. For 7 = 12, we have A4(n = 265) = 4.38238, M(n = 512) = 
4.38589. 

As was remarked in [2], for the value of 7 used in [3], the PCSM mass gap, m, is very 
small, compared to the coupling e, and becomes smaller as 7 is increased. This means that 
it is very difficult to observe the confined phase, since string breaking occurs very easily. 
Since we have shown that there exist finite-mass gluons in the lattice theory at finite volume, 
it is reasonable to observe a Higgs-like potential between two sources at large separations, 
which explains the smooth crossover of [3]. 
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We have shown that at infinite N, the gluon mass has qualitatively similar behavior to 
the one observed numerically in [3] for N = 2. Now the very important question remains, 
what can we actually say analytically about the finite -N case (particularly N = 2)1 

Reproducing a nonperturbative result like the correlation functions (2.4), (3.2) for 
general N is a very difficult task that we are unable to achieve at this point. There are 
several difficulties to overcome. 

First of all, the form factors of the PCSM renormalized field at finite N are not known, 
and it is significantly harder to calculate them. In fact, the simplicity of the S-matrix at 
large N is what made the calculation of form factors possible. A notable exeption is the 
N = 2 case, where the PCSM is equivalent to an 0(4)-symmetric nonlinear sigma model, 
using the fact that SU(2) x SU(2) ~ 0(4). Some of the few-particle form factors of the 
0(4) model have been calculated in [14], However, we need an explicit expression of all the 
form factors to study the ultraviolet regime of the correlation function, so the results of [14] 
are not enough to complete the N = 2 computation. 

Once all the form factors are known there is still the second problem of extracting the 
ultraviolet information from the correlation function by analyzing the spectral sum (2.1) 
as was done in [5] for large N. This was a very nontrivial task that was possible only 
because of the simple structure of the large -N form factors, where 0(1/N) corrections were 
discarded at various points. It will very likely take a lot more time and effort to extract a 
simple result like (2.4) from the finite-A r form factors, even if the form factors were known. 

Finally, if we want to study the finite -N correlation function at finite volume/temperature, 
the simple Leclair-Mussardo formula is no longer valid, and it is not clear what is the way 
to proceed with this calculation. As was shown in [6], the Leclair-Mussardo formula is only 
expected to be valid at infinite N, because of two important properties: the expansion is 
well defined (there are no poles in the real line of rapidities), and the thermal distributions 
arising from the thermodynamic Bethe ansatz are trivial (like that of an ideal gas). There 
is no reason to believe these properties persist at finite N. It is likely that a more careful 
regularization scheme is needed, like the one proposed in [15]. 

Having established the difficulty of extending our results to general N. one can consider 
the possibility of studying large, but not infinite N. by finding small corrections to our result 
in powers of 1/N. In principle this could be possible, even though it might require a very 
difficult and careful calculation, since there are several steps of the computation where 
O (1/A^) corrections have been ignored. To compute these small corrections, one needs to 
keep track of the terms that have been ignored at the levels of the S-matrix, computation 
of form factors, TBA, and correlation functions. 

At the S-matrix level, it is very simple to find 1/N corrections, since and exact ex¬ 
pression is known for all N [8]. However, the most important simplifying property of the 
infinite -N limit is that, as was discussed in [6], the scattering becomes effectively diagonal. 

In this case this means that most particles don’t interact with each other, unless there is 
a color-index contraction between them. Small 1 /N corrections already present a difficulty 
when computing form factors, since it means that all particles interact with each other. 
There is no fundamental reason why this computation should not be possible, but there 
will be many more contributions to the form factors to take into account from new particle 
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interactions that were not present at infinite N. 

Onr TBA and Leclair-Mussardo correlation-function computations were significantly 
simplified by relying on the fact that the scattering was diagonal. TBA computations are 
much more difficult in a non-diagonal theory, since one needs to take into account the 
thermodynamic contributions of pseudo particles, such as magnons and their bound states 
(also called "strings"). A full non-diagonal TBA computation for the PCSM was proposed 
in [16], but the resulting functional equations have not yet been solved for large values of 
N. The original proposal for computing correlations functions by Leclair and Mussardo 
was also only applicable to diagonal scattering theories. A proposal for generalizing the 
computation of these correlation functions to non-diagonal theories was made in Ref. [17]. 

There is in principle not a physical reason why it should not be possible to compute 1 /N 
corrections. But as we have discussed, even a small correction would completely change the 
techniques we have to use, since the PCSM is a diagonal scattering theory only at precisely 
N = oo. 

Despite all these negative arguments, there are still some statements we can make about 
the Higgs phase at finite N. There are, however, two disadvantages: we are forced to rely 
only on perturbative calculations, and we can only study the infinite volume case. 

A simple perturbative analysis tells us how the general -N two-point function diverges 
at x = 0. This simple result can be found in the first reference of [11], and is also discussed 
in [5]. From leading-order perturbation theory one can find the time-ordered two-point 
function, G(x, A), and the coupling go(A) satisfy the renormalization group equations: 


dlnG(x, A) 
d In A 


7(9o) = 7iffo + ■ ■ ■ 


dg 2 0 (A) 

Sin A 


P(do) 


~Pi9o + ‘ j 


(4.3) 


with the coefficients 71 = ( N 2 — 1)/(2itN 2 ) and f3\ = 1/47T. Integrating (4.3), one finds 


G(0,A) 


C 


In [ — 

m 


71/ft 


+ ' 


(4.4) 


where C is some undetermined constant. Equation (4.4) is completely consistent with the 
nonperturbative (2.4), since liniAr^oo 7 i// 3 i = 2, which is the main result of Ref.[5]. One 
advantage of the nonperturbative calculation is that one can find the constant C = C 2 , 
while for finite N this is unknown from this calculation. 

This very simple perturbative result allows us to write an expression for any N for the 
gluon mass at infinite volume, up to the constant C. Combining (4.4) with (4.1), we can 
write for N = 2 (using the notation of Eq. (4.2), for n — 00 ) 


M( 7 ,/?) 



In 



1 

2 


(4.5) 


Equation (4.5) is not good enough to compare with the results from [3] because as we 
argued, their computations are done in very small volumes. At small volumes we expect 
there might appear some volume-dependent contributions like those of (4.2) that cannot be 
ignored. 
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Even though the result (4.5) is only valid for infinite volume, it can still show the basic 
conclusion of this paper. That is, introducing a lattice regularization makes the gluon mass 
finite, and even small, depending on the values of the couplings. For example, we show 
some values that span the range used in [3], .Ad (2,120) = 0.843767\/C', .Ad (10,120) = 
6.35361^, M( 2, 200) = 0.653579\/C, M(8, 200) = 3.71684v / C. 
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